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Introduction:  Prostate  cancer  (PCa)  is  the  most  common  cancer  in  men  in  several 
countries,  with  the  American  Cancer  Society  (ACS)  estimating  241,740  new  cases  of  PCa 
to  be  diagnosed  and  deaths  of  28,170  men  of  PCa  (1).  Due  to  its  prevalence  in  the  male 
population  as  well  as  its  unpredictable  clinical  course,  early  detection  and  diagnosis  have 
become  a  priority  for  many  health  care  professionals.  Another  method  for  staging 
prostate  cancer  is  through  imaging  techniques  including  ultrasound,  computed 
tomography  (CT),  and  magnetic  resonance  imaging  (MRI)  with  or  without  the  help  of 
dynamic  contrast  enhancement  modeling  (DCE-MRI),  diffusion  weighted  imaging 
(DWI),  and  magnetic  resonance  spectroscopy  (MRS)  (2-5).  MRS  is  a  powerful  tool  for 
exploring  the  cellular  chemistry  of  human  tissues  (3,5,6-11).  There  is  a  growing  body  of 
evidence  that  1H  MRS  may  contribute  to  the  clinical  evaluation  of  prostate  cancer  and 
also  for  evaluating  the  metabolic  alterations  due  to  therapy.  There  have  been  no  reports 
on  combining  two  spectral  dimensions  with  two-dimensional  (2D)  or  three  dimensional 
(3D)  spatial  encoding  applicable  to  prostate  cancer.  Acceleration  of  magnetic  resonance 
spectroscopic  imaging  (MRSI)  has  been  demonstrated  using  echo-planar  imaging 
techniques  (12-13).  Recently,  Schulte  et  al.  have  successfully  developed  an  algorithm 
called  prior-knowledge  fitting  (ProFit)  to  quantify  metabolite  concentrations  using  the 
JPRESS  spectra  recorded  using  a  Philips  3T  MRI  scanner  (14).  It  was  demonstrated  that 
metabolite  quantitation  of  JPRESS  spectra  with  ProFit  was  accurate,  robust  and  yielding 
generally  consistent  results,  both  in  vivo  and  in  vitro.  Their  results  suggest  that  the 
number  of  quantifiable  prostate  metabolites  can  be  increased  from  3-4  with  ID 
PRES S/LC -Model  to  more  than  10  with  JPRESS/ProFit  (15-16). 

Body: 

i)  Proposed  Task  1  ( Months  1-6):  To  implement  a  multi-voxel  based  extension  of  the 
JPRESS  sequence,  in  which  two  spectral  encodings  will  be  combined  with  two  spatial 
encodings  using  the  new  Siemens  VB17a  platform.  This  four-dimensional  (4D)  data 
acquisition  scheme  will  be  accomplished  utilizing  the  EPI  approach  that  is  commonly 
used  for  spatial  encoding  in  MRI. 

Completed  and  Reported  in  the  1st  year  Annual  Report. 

ii)  Proposed  Task  2:  To  evaluate  the  EPI-based  JPRESS  using  a  prostate  phantom 
containing  several  metabolites  which  have  been  reported  in  prostate  tissues,  and  to 
optimize  the  EP-JJRES1  sequence  and  other  acquisition  parameters  using  the  phantom 

(Months  6-12 ). 

Completed  and  Reported  in  the  1st  year  Annual  Report. 

iv)  Proposed  Task  3:  To  record  the  4D  EP-JRES1  spectra  in  the  peripheral,  central  and 
transition  zones  of  healthy  prostates.  (Months  6-18). 

Five  more  healthy  males  have  been  investigated. 

Proposed  Task  4:  To  develop,  evaluate  and  optimize  the  prior-knowledge  basis  set 
spectra  using  the  GAMMA- simulation  and  prostate  phantom  solutions  as  prior 
knowledge  for  the  multi-voxel  based  JPRESS  spectra  recorded  using  the  3T  MRI  scanner 

(Months  6-12). 
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Accomplished  during  September  29, 2011-October  28  2012: 


Proposed  Task  5:  To  record  the  multi-voxel-based  2D  spectra  in  the  peripheral,  central 
and  transition  zones  of  patients  with  BPH  and  malignant  prostate  cancer.  The  prostate 
metabolite  concentrations  calculated  using  the  ProFit  algorithm  prepared  for  the  multi¬ 
voxel  data  will  be  compared  with  LC-Model  processed  ID  spectral  based  MRSI  data 

(Months  18-36). 


Accomplished  during  October  29,  2013-October  28  2014:  Six  malignant  prostate 
cancer  patients  and  five  healthy  males  were  investigated  during  the  current  year.  The  4D 
EP-JRESI  data  were  recorded  using  the  following  parameters:  TR/TE=1.5s/30ms,  2 
averages,  512  t2,  oversampled  32kx,  64  increments  along  the  indirect  spectral  (ti)  and  16 
spatial  ky  dimensions;  the  endorectal  "receive"  coil.  The  oversampled  32kx  were 
processed  to  normal  16kx  points.  First,  the  4D  NUS  EP-JRESI  data  was  acquired  in  a  71 
year  old  PCa  patient;  Extracted  1ml  spectra  from  3  different  locations  after  the  CS 
reconstruction  are  shown  in  Figure  1  (ii,  iii  and  iv)  and  the  middle  insert  (i)  shows  the 
multi-voxel  grids  overlaid  on  an  T2W  axial  MRI  slice. 


3.8  3.6  3.4  3.2  3.0  2.8  2.6  2.4  2.2  2.0 
F2  (ppm) 


3.8  3.6  3.4  3.2  3.0  2.0  2.6  2.4  2.2  2.0 
F2  (ppm) 


ii)  Cancer 


iii)  Cancer 


F,=0(Hz) 


Figure  1.  i)  Axial  T2W  slice  showing  the  multi-voxel  EP-JRESI  grids;  extracted  2D 
JPRESS  spectra  of  malignant  voxels  (ii)  and  iii),  and  the  healthy  voxel  in  the  peripheral 
zone  (iv). 
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We  were  able  to  reproduce  the 
above  mentioned  findings  in  more 
patients  with  known  GS,  and  also, 
in  remaining  four  patients  whose 
histology  findings  were  not 
known.  2D  peak  (cross  and 
diagonal)  volumes  were  defined 
by  the  operator  in  the  frequency 
domain  2D  JPRESS  data  for  this 
pilot  analysis.  Two  different  ratios 
were  calculated  using  the  4D 
NUS  EP-JRESI  data  acquired  in 

Glx  Glx  Cit  Cit  1 

(3+3/3+4)  (4+3/4+4)  <3+3/3+4)  (4+3/4+4)  14  PCa  patients: 

(Glx/(Ch+Cr+Spm)  and 

Cit/(Ch+Cr+Spm),  one  set  from  the  affected  lesions  and  another  set  from  the  contralateral 
healthy  sides.  As  evident  in  Fig.2,  additional  metabolite  (Glx)  changes  can  be  detected 
using  the  4D  NUS  EP-JRESI  data.  Both  ratios  showed  significant  changes  (p<0.05); 
however,  this  outcome  from  a  small  group  of  PCa  patients  needs  to  be  treated  with 
caution  showing  a  necessity  for  further  research  using  larger  patient  cohorts. 

Key  Research  Accomplishments 

•  Continued  evaluation  of  the  4D  EP-JRESI  sequence  in  6  malignant  prostate 
cancer  patients  and  five  healthy  male  subjects. 

•  We  have  made  further  progress  on  compressed  sensing  reconstruction  of  the 
non-uniformly  undersampled  4D  EP-JRESI  sequence  using  different 
reconstruction  methods  such  as  maximum  entropy,  total  variation  (TV),  etc  .  This 
will  facilitate  shortening  the  endorectal  spectral  acquisition  and  reducing  the 
patient  inconvenience  during  the  scan.  A  manuscript  was  submitted  to  NMR  in 
Biomedicine  and  a  major  revision  is  in  progress. 

•  Using  the  preliminary  results  obtained  using  this  IDEA  grant,  we  had  reported  our 
NIH  R01  grant  application  submission  entitled  “Fast  J-resolved  Prostate  MR 
Spectroscopic  Imaging  and  Non-linear  Reconstruction”  was  resubmitted  on  July 
5,  2014  using  the  preliminary  results  from  this  IDEA  grant.  The  percentile  score 
was  37%  and  we  plan  to  resubmit  this  grant  to  NIH. 

Reportable  Outcomes: 

A.  Peer-reviewed  Publications: 

1.  Nagarajan  R,  Iqbal  Z,  Burns  B,  et  al.  Accelerated  Echo  Planar  J-Resolved 
Spectroscopic  Imaging  in  Prostate  Cancer:  Nonlinear  Reconstruction  Using  Total 
Variation  And  Maximum  Entropy.  NMR  in  Biomedicine  2014  [being  revised]. 


O 

+ 


r  i.5 


Figure  2.  Glx  and  Cit  ratios 
from  7  (GS3+4,  GS4+3)  and 
7  ( GS4+3 ,  GS4+4)  PCa 
patients. 


B.  Presentations:  During  the  3rd  year,  three  abstracts  were  presented  at  the  23rd 
International  Society  of  Magnetic  Resonance  in  Medicine  (ISMRM),  Milan,  Italy,  May 
10-16,  2014:  the  first  abstract  entitled  “Accelerated  Echo  Planar  J-Resolved 
Spectroscopic  Imaging  in  Prostate  Cancer:  Nonlinear  Reconstruction  Using  Total 
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Variation  And  Maximum  Entropy”  (oral  presentation),  #962;  the  2nd  abstract  entitled 
“Maximum  Entropy  Reconstructed  Echo  Planar  J-Resolved  Spectroscopic  Imaging  and 
Diffusion  Weighted  Imaging  in  Prostate  Cancer”  (eposter),  #  4108;  the  3rd  abstract 
entitled  “Compressed  Sensing  based  Echo-planar  3D  MRSI  Using  Short  Echo  Time:  A 
Pilot  Evaluation  Using  a  Prostate  Phantom”  (poster),  #  1529  summarizing  the 
implementation  of  a  novel  3D  EPSI  sequence  and  evaluation  of  it  using  a  prostate 
phantom. 

C.  Books:  None  on  Prostate  Cancer  Research  based. 

Conclusions:  After  successfully  implementing  the  4D  EP-JRESI  scanning  protocol  on 
the  3T  MRI  scanner  and  testing  it  in  healthy  controls  and  patients  during  the  earlier  2 
years,  11  more  subjects  have  been  investigated  in  the  3rd  year.  We  will  continue  to 
recruit  10  malignant  and  5  benign  prostate  cancer  patients,  and  5  healthy  male  subjects 
during  the  next  year  (no-cost  time  extension). 
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ACCELERATED  ECHO  PLANAR  J-RESOLVED  SPECTROSCOPIC  IMAGING  IN  PROSTATE  CANCER:  NONLINEAR 
RECONSTRUCTION  USING  TOTAL  VARIATION  AND  MAXIMUM  ENTROPY 

Rajakumar  Nagarajan1,  Zohaib  Iqbal1,  Brian  Burns1,  Neil  Wilson1,  Manoj  K  Sarnia1,  Daniel  A  Margolis1,  Robert  E  Reiter2,  Steven  S  Raman1,  and  M.Albert  Thomas1 
Radiological  Sciences,  University  of  California  Los  Angeles,  LOS  ANGELES,  CA,  United  States,  2  Urology,  University  of  California  Los  Angeles,  LOS  ANGELES,  CA, 

United  States 


Target  audience:  MR  pulse  sequence  developers,  scientists/clinicians  interested  in  accelerated  imaging  and  non-linear  reconstruction  in  prostate  imaging 


Introduction:  Prostate  cancer  is  the  second  leading  cause  of  cancer  death  in  men  next  to  lung  cancer.  The  American  Cancer  Society's  estimates  238,590  new 

cases  and  29,720  deaths  from  prostate  cancer  in  2013  (1).  One  in  6  men  will  get  prostate  cancer  during  his  lifetime  and  one  in  36  will  die  of  this  disease.  Current 
limitations  of  single  voxel  and  multivoxel  prostate  spectroscopic  imaging  are  due  to  the  overlap  of  metabolite  resonances,  quantifying  few  metabolites  only  (citrate 
(Cit),  choline  (Ch),  creatine  (Cr)  and  spermine  (Spm))  and  long  echo  times.  A  single- voxel  based  two-dimensional  (2D)  J  resolved  spectroscopic  sequence  (JPRESS) 
has  been  implemented  and  demonstrated  in  prostate  cancer,  and  showed  improved  spectral  dispersion  due  to  the  added  spectral  dimension  (2,3).  Compressed  sensing 
was  recently  introduced  as  a  powerful  method  to  accelerate  MRI  by  exploiting  the  sparsity  of  the  images  in  a  known  transform  domain  to  reconstruct  non-uniformly 
undersampled  (NUS)  k-space  data  (4).  By  combining  EPSI  with  JPRESS,  2D  spectra  can  be  recorded  in  multiple  locations  in  prostate  using  four  dimensional  (4D) 
Echo-Planar  J-Resolved  Spectroscopic  Imaging  (EP-JRESI)  combining  2  spectral  with  2  spatial  dimensions.  A  pilot  feasibility  was  demonstrated  recently  to  map 
metabolites  in  the  human  prostate  (5).  Maximum  entropy  (MaxEnt)  and  total  variation  (TV)  image  reconstruction  have  been  used  to  reconstruct  NUS  indirect  spectral 
and  spatial  dimensions  (6).  The  goal  of  the  present  work  was  to  compare  and  correlate  the  significant  metabolites  of  the  accelerated  4D  EP-JRESI  data  acquired  in 
prostate  cancer  and  reconstructed  with  MaxEnt  and  TV. 

Materials  and  Methods:  Twenty  two  prostate  cancer  patients  with  a  mean  age  of  63.8  years  (range:  46-79  years)  were  investigated  in  this  study.  Prostate- 

specific  antigen  (PSA)  varied  from  0.7  to  22.8  ng/mL  (mean  of  6.23  ng/mL)  and  the  mean  delay  between  biopsy  and  MR  investigation  was  8  weeks.  This  study  was 
approved  by  the  Institutional  Review  Board,  and  informed  consent  was  obtained  from  each  patient.  A  Siemens  3T  MRI  Scanner  with  high-performance  gradients  (Trio- 
Tim,  Siemens  Medical  Solutions,  Erlangen,  Germany)  was  used  in  this  investigation.  A  4D  NUS  EP-JRESI  sequence  was  implemented  on  the  3T  MRI  scanner  and  the 
volume  of  interest  (VOI)  was  localized  using  three  slice-selective  radio-frequency  (RF)  pulses  (90°-180°-1800).  The  total  time  for  obtaining  a  fully  sampled  4D  EP- 
JRESI  data  (TR  of  1.5s,  16ky*16kx,  64-100ti,  512t2)  can  be  more  than  25  minutes.  The  MaxEnt  and  TV  reconstructions  enable  a  stable  and  accurate  reconstruction 
from  the  NUS  based  EP-JRESI  data  (25%  of  ti  and  ky  increments)  with  a  reduced  total  time.  The  parameters  for  the  EP-JRESI  were:  TR/TE/Avg  =  1.5s/30ms/l-2,  16 
phase  encoding  steps,  512  complex  points  with  an  F2  bandwidth  of  1190Hz.  For  the  second  dimension  (Fi),  64  increments  with  bandwidths  of  1000Hz  were  used.  A  4X 
NUS  was  imposed  along  the  plane  containing  incremented  spectral  and  spatial  dimensions  (ti  and  ky).  The  individual  voxel  volume  in  human  prostate  was  1ml.  Two 
sets  of  data  were  collected,  one  with  water  suppression  (WS)  with  a  total  scan  time  of  6-12  minutes  and  second  with  non- water  suppression  (NWS)  using  1  average 
only.  The  NUS  data  was  reconstructed  by  MaxEnt  and  TV  separately.  A  modified  Split-Bregman  algorithm  (6)  solves  the  unconstrained  TV  optimization  problem  as 
shown  in  Eqn.(l)  below: 

min||Vm||i  +X\\  Fum-y||2  (1) 

m 


where  V  is  the  gradient  operator,  m  is  the  reconstructed  data,  ||x||i  is  the  f-norm,  X  is  a  regularization  parameter,  Fu  is  the  under  sampled  Fourier  transform,  and  y  is  the 
under-sampled  data.  The  above  equation  removes  the  incoherent  artifacts  due  to  NUS  by  minimizing  the  total  variation  (TV)  while  maintaining  consistency  with  the 
sampled  measurements.  MaxEnt  is  a  constrained  convex  optimization  algorithm  that  uses  a  variant  of  the  conjugate  gradient  method  to  iteratively  solve  the  inverse 
problem  (6-8): 

maximize  Si/2  (f)  s.t.  ||F_1  Kf-D||2  <  o  (2) 

where/is  the  estimated  fully-sampled  spectrum  at  each  iteration,  F1  is  the  inverse  Fourier  transform,  K  is  the  NUS  matrix,  D  is  the  measured  time-domain  data,  o  is  the 
noise  standard  deviation,  and  S1/2(f)  is  the  spin-1^  entropy  of  the  estimated  spectrum  (6-8). 

Results  and  Discussion:  All  NUS  based  4D  EP-JRESI  data  were  processed  using  TV  and  MaxEnt  reconstruction  with  custom  processing  software.  We 

observed  decreased  levels  of  Cit,  Spm  and  ml  with  respect  (Ch+Cr)  in  cancer  locations  compared  to  non-cancer  locations  in  both  TV  and  MaxEnt  reconstruction.  Also 
we  have  compared  and  correlated  the  TV  with  MaxEnt  reconstruction  methods  of  a  few  metabolites  which  showed  significant  decrease  as  well  as  decreasing  trend  in 
the  cancer  locations  compared  to  non-cancer  locations.  Figure  1  shows  the  correlation  of  maxEnt  and  TV  of  Cit/(Ch+Cr)  and  Spm/(Ch+Cr)  in  cancer  and  non-cancer 
location.  Similarly  a  positive  correlation  was  found  in  the  following  metabolites  in  cancer  locations:  Cit/(Ch+Cr)  (R2=0.8708),  mE(Ch+Cr)  (R2=0.972)  and 
Spm/(Ch+Cr)  (R2=0.7348).  Also  we  observed  positive  correlation  of  Spm  and  ml  in  non-cancer  location.  Table.  1  shows  the  metabolites  ratio  of  maxEnt  and  TV  in 
cancer  and  non-cancer  locations.  There  was  an  overlap  between  cancer  and  non-cancer  locations,  possibly  due  to  the  low  signal-to-noise  ratio  of  the  Cr  peak.  In 
addition,  Ch  and  Cr  were  not  always  well  resolved  especially  in  cancer  locations,  adding  to  the  uncertainty.  The  advantage  of  NUS  based  4D  EP-JRESI  sequence  will 
record  short  TE-based  spectra  from  multiple  regions  of  a  prostate  and  will  show  more  metabolites,  such  as  glutamate  (Glu),  glutamine  (Gin),  myo-inositol  (ml),  taurine 
(Tau),  lactate  (Lac),  and  choline  groups,  in  addition  to  the  normally  detected  Cit,  Cr,  Ch  and  Spm. 


Spm/(Ch+Cr) 


maxEnt 

Figure  1.  Correlation  of  MaxEnt-processed  data  with  TV  of  Cit  and  Spm  ratios 


Metabolites 

Cancer 

Non-Cancer 

wrt  (Ch+Cr) 

(Mean  ±  SD) 

(Mean  ±  SD) 

maxEnt 

Cit 

0.378+0.213 

1.115±0.442 

Spm 

0.498+0.323 

0.633+0.560 

ml 

0.189±0.119 

0.232±0.146 

TV 

Cit 

0.413±0.252 

1.108+0.449 

Spm 

0.505+0.319 

0.918±0.817 

ml 

0.214+0.154 

0.276±0.215 

Table.l.  Metabolites  ratio  of  maxEnt  and  TV 


Conclusion:  We  were  able  to  detect  and  quantity  metabolites  using  NUS  based  4D  EP-JRESI  data  acquired  in  clinically  acceptable  time  (<12  minutes).  We 

have  shown  that  it  possible  to  undersample  the  4D  EP-JRESI  sequence  with  an  acceleration  factor  of  4X  and  the  data  can  be  reliably  reconstructed  using  TV  and 
MaxEnt  methods.  Both  non-linear  reconstruction  methods  provide  comparable  results. 
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Introduction:  In  the  United  States,  90%  of  men  with  prostate  cancer  (PCa)  are  older  than  60  years,  diagnosed  by  early  detection  with  the  serum  prostate-specific 
antigen  (PSA)  blood  test,  and  have  disease  believed  to  be  confined  to  the  prostate  gland  (1).  In  vivo  proton  MRS  provides  biochemical  information  on  the  prostate 
tissues  by  providing  relative  concentrations  of  citrate  (Cit),  creatine  (Cr),  choline  (Ch)  and  polyamines  that  are  used  to  detect  and  diagnose  PCa  (2).  The  challenging 
task  in  ID  MRS  spectra  of  PCa  is  the  quantification  of  individual  metabolites,  due  to  spectral  overlap  of  metabolites  because  of  limited  spectral  dispersion  at  clinically 
used  magnetic  field  strengths  of  less  than  3T.  Conventional  phase  encoded  MRSI  is  relatively  inefficient  and  time-consuming  because  it  involves  a  large  number  of 
phase  encodings.  The  applicability  of  compressed  sensing  (CS)  (3)  to  four  dimensional  (4D)  echo  planar  based  J-Resolved  spectroscopic  imaging  (EP-JRESI)  in  order 
to  accelerate  data  acquisition  using  nonuniform  under  sampling  (NUS)  (4).  This  could  be  significantly  reduced  without  sacrificing  the  spectral  quality  while  maintaining 
high  spatial  resolution  as  well  as  the  4D  EP-JRESI  detect  more  metabolites  (ml,  glutamate,  glutamine,  spermine  and  lactate)  compared  to  long  TE  conventional  method. 
In  this  study  maximum  entropy  (MaxEnt)  method  was  used  to  reconstruct  the  4D  EP-JRESI  data.  Diffusion  weighted  imaging  (DWI)  has  the  ability  to  qualitatively  and 
quantitatively  represent  the  diffusion  of  water  molecules  by  the  apparent  diffusion  coefficient  (ADC),  which  indirectly  reflects  tissue  cellularity  (5).  The  major  goal  of 
the  study  is  to  demonstrate  the  feasibility  of  correlating  MaxEnt  reconstructed  4D  EP-JRESI  data  with  DWI-MRI  in  PCa  patients  in  3T. 

Materials  and  Methods:  4D  EP-JRESI  and  DWI )  were  recorded  in  eighteen  PCa  patients  ranging  in  age  from  46-73  years  (mean,  65  years  using  Siemens  3T  MRI 
Scanner  (Siemens  Medical  Systems,  Germany).  This  study  was  approved  by  the  Institutional  Review  Board,  and  informed  consent  was  obtained  from  each  patient.  The 
mean  delay  between  biopsy  and  MRI  was  8  weeks.  A  4D  NUS  EP-JRESI  sequence  was  implemented  on  a  Siemens  3T  Trio-Tim  scanner  and  the  volume  of  interest 
(VOI)  was  localized  using  three  slice-selective  radio-frequency  (RF)  pulses  (90°-180°-1800).  The  CS-based  reconstruction  will  enable  a  stable  and  accurate 
reconstruction  from  the  NUS  based  EP-JRESI  data  (25%  of  fi  and  ky  increments)  with  a  reduced  total  time.  The  parameters  for  the  EPJRESI  were:  TR/TE/Avg  = 
1.5s/30ms/l-2,  16  phase  encoding  steps,  512  complex  points  with  an  F2  bandwidth  of  1190Hz.  For  the  second  dimension  (FI),  64  increments  with  bandwidths  of 
1000Hz  were  used.  A  25%  NUS  was  imposed  along  the  incremented  spectral  and  spatial  dimensions.  The  individual  voxel  volume  in  human  prostate  was  1ml.  MaxEnt 
is  a  constrained  convex  optimization  algorithm  that  uses  a  variant  of  the  conjugate  gradient  method  to  iteratively  solve  the  inverse  problem  (6): 


maximize  Si/2  (f)  s.t.  HF1  Kf-D||2  <  o  (1) 

where  f  is  the  estimated  fully-sampled  spectrum  at  each  iteration,  F-l  is  the  inverse  Fourier  transform,  K  is  the  NUS  matrix,  D  is  the  time-domain  measured  data,  o  is 
the  noise  standard  deviation,  and  Sl/2(f)  is  the  spin-1/!  entropy  of  the  estimated  spectrum  (7).  In  addition,  DWI-MRI  was  acquired  with  the  following  parameters: 
(TR/TE  3200/60  ms,  bandwidth  1250  Hz  in  the  EPI  frequency  direction)  with  b  values  of  0,  100,  400  and  800  s/mm2. 

Receiver  operating  characteristic  (ROC)  curve  analyses  based  on  logistic  regression  models  were  performed  to  predict  the  MRSI  and  ADC  values.  P  values  <  0.05  were 
considered  statistically  significant. 

Results:  Figure  1  show  the  ratios  of  myo-inositol  (ml),  spermine  (Spm)  and  Cit  with  respect  to  Cr  as  well  as  (Ch+Cr).  Fig.2  shows  the  ADC  values  between  cancer  and 
non-cancer  locations  of  PCa  patients.  An  investigation  of  the  measurement  of  EP-JRESI  with  DWI  data  showed  interesting  results  in  this  study.  Since  Ch  and  Cr  are 
not  always  well  resolved,  adding  to  the  uncertainty,  we  have  calculated  two  different  ratios  (Cr,  (Ch+Cr)).  Table  1  shows  the  receiver  operating  characteristic  (ROC) 
curve  analyses  based  on  logistic  regression  models  were  performed  to  predict  the  EP-JRESI  and  ADC  values.  Combining  both  EP-JRESI  and  ADC  gives  100% 
accuracy  in  this  selected  group  of  patients.  We  have  detected  myo-inositol,  glutamine,  glutamate,  spermine  in  addition  to  Cit,  Cr  and  Ch.  The  NUS  based  4D  EP-JRESI 
data  were  processed  using  MATLAB  based  post  processing.  We  have  found  significant  (p<  0.05)  decrease  of  Cit,  spm,  with  respect  to  (Ch+Cr)  in  cancer  locations 
compared  to  non-cancer  locations.  Also  decreased  trend  of  ml  was  observed  with  respect  to  (Ch+Cr)  ratio.  Similar  trends  followed  in  Cr  ratio  except  ml.  However, 
because  of  the  small  number  of  samples  in  each  group  the  differences  in  ml/Cr  ,  mI/(Ch+Cr)  did  not  reach  statistical  significance.  The  ROI  based  mean  ADC  values  of 
cancer  and  non-cancer  locations  in  the  peripheral  zones  were:  (0.985+0.164)  and  (1.453+0.168).  The  statistical  significance  (p  <0.05)  was  observed  between  cancer  and 
non-cancer  locations. 


Metabolites 


Fig  1.  Metabolites  ratios  with  respect  to  (Ch+Cr) 
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Fig  2.  ADC  values  from  Cancer  and  non-cancer 


Parameter  (%) 

MRSI 

ADC 

Sensitivity 

86.7 

95.2 

Specificity 

95.5 

95.0 

PPV 

91.3 

94.3 

NPV 

92.8 

93.0 

Accuracy 

91.5 

93.1 

AUC 

97.0 

97.4 

Table.l.  Findings  from  EP-JRESI  and  DWI 


Discussion:  Normal  prostate  epithelial  cells  produce  and  accumulate  a  large  amount  of  Cit  which  is  secreted  as  a  major  component  of  the  prostatic  fluid.  Compared  to 
normal  tissue,  decreased  levels  of  Cit  are  previously  observed  in  PCa  tissue  by  ex  vivo  MRS  (8). why  ex  vivo  only  Decrease  in  polyamines  is  associated  with  prostate 
cancer  (9).  Additionally,  the  very  low  putrescine  concentration  in  our  study  confirms  that  the  polyamine  peak  predominantly  consists  of  Spm.  Spm  cannot  be  fully 
separated  from  the  Ch  peak  using  ID  spectral  based  MRSI,  but  with  EP-JRESI  we  were  able  to  detect  and  quantify  nicely.  The  increase  of  ml/Cr  was  discussed  in  PCa 
earlier  which  agreed  with  our  result  (10)  and  decrease  of  mI/(Ch+Cr)  may  be  the  overlap  of  Ch  and  Cr  which  needs  to  be  investigated  in  the  future.  The  lower  mean 
ADC  values  in  the  more  aggressive  tumors  in  our  study  may  be  due  to  higher  cellular  density  in  poorly  differentiated  tumors,  resulting  in  more  restricted  movement  of 
water  protons  (11). 

Conclusion:  It  is  evident  from  our  preliminary  investigation  using  a  3T  MRI  scanner  that  NUS  based  4D  EP-JRESI  with  MaxEnt  reconstruction  facilitates 
unambiguous  detection  of  less  T2- weighted  metabolites  compared  to  conventional  ID  MRSI.  DWI  is  potentially  useful  in  PCa  detection  and  localization.  Our  pilot 
findings  demonstrate  a  possible  clinical  application  of  the  4D  NUS  EP-JRESI  sequence  in  combination  with  DWI  which  may  be  useful  in  localizing  PCa  for  diagnosis, 
monitoring,  biopsy,  treatment  planning  and  staging  as  well  as  potentially  for  response  assessment. 
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Introduction:  MRSI  of  the  prostate  is  typically  performed  with  a  combination  of  point-resolved  spectroscopy  (PRESS)  (1)  volume  localization  and  three  dimensional 
(3D)  MRSI  (2)  rather  than  the  traditional  single-voxel  or  two-dimensional  (2D)  MRSI  technique  used  for  many  years  in  brain  imaging.  Three-dimensional  MRSI 
requires  phase  encoding  in  three  dimensions,  conventionally  known  as  frequency,  phase,  and  slice.  Acquisition  time  and  coverage  of  the  prostate  are  the  main 
considerations  in  choosing  the  matrix  dimensions.  The  MRSI  data  can  be  acquired  with  higher  spatial  resolution  leading  to  long  total  acquisition  time.  Conventional 
prostate  MRSI  study  with  average  weighted  encoding  uses  long  echo  time  (TE)  with  short  repetition  time  (TR)  allowing  the  observation  of  a  reduced  number  of 
metabolites  without  accurate  quantitation.  Long  TEs  are  used  for  the  3D  MRSI  due  to  the  addition  of  MEGA  radio-frequency  pulses  for  both  water  and  lipid 
suppression.  In  echo-planar  spectroscopic  imaging  (EPSI)  (3,4),  an  alternating  gradient,  which  simultaneously  encodes  space  and  chemical  shift  (time),  is  applied  along 
one  of  the  spatial  (readout)  directions,  thereby  reducing  phase  encoding  dimensions  from  three  to  two  thereby  increasing  speed.  MRI  is  well  suited  for  compressed 
sensing  (CS),  and  there  are  significant  benefits  in  imaging  speed  and  reduced  costs,  thereby  improving  patient  care  (5,6).  A  major  challenge  in  designing  CS  data 
acquisition  methods  for  MRI  is  in  implementing  non-uniform  under  sampling  (NUS)  densities  that  result  in  incoherent  aliasing  while  providing  data  sparsity  in  a 
transform  domain,  such  as  wavelets,  curvelets,  etc.  The  goal  of  the  present  study  was  to  implement  a  novel  NUS  based  3D  EPSI  on  a  3T  MRI/MRS  scanner  and  to 
evaluate  the  performance  using  a  prostate  phantom. 

Materials  and  Methods:  A  schematic  diagram  of  NUS  based  3D  EPSI  sequence,  compiled  and  implemented  on  the  Siemens  3T  MRS  scanner,  is  shown  in  Figure  1 
where  three  spatial  encodings  and  one  spectral  dimension  were  used.  In  the  3D  EPSI  sequence,  the  EPSI  readout  simultaneously  acquires  one  spatially  encoded 
dimension  (kx)  and  one  spectral  dimension  (t),  leaving  the  remaining  the  two  spatial  dimensions  (ky  and  kz,  respectively)  incremented  sequentially.  We  propose  the  use 
of  NUS  in  the  remaining  ky-kz  plane,  using  CS  to  reconstruct  the  equivalent  missing  data  to  a  fully  sampled  3D  EPSI  acquisition. 

To  determine  the  feasibility  as  well  as  evaluate  the  performance  of  the  CS  reconstruction,  numerous  prospective  2X  undersampling  and  reconstructions 
were  performed  on  a  prostate  3D  EPSI  of  nine  phantom  dataset  and  two  healthy  volunteers  using  external  phased  array  coil.  A  500mL  prostate  phantom  was  prepared 
containing  the  following  metabolites  at  physiological  concentrations  as  reported  in  healthy  human  prostate  (7):  Cit,  50  mM,  Cr,  5  mM,  Ch,  1  mM,  Spm,  6  mM,  ml,  10 
mM,  PCh,  2  mM,  Tau,  3  mM,  Glu,  4  mM,  Gin,  2.5mM  and  si,  0.8  mM.  The  dataset  was  localized  with  a  field  of  view  (FOV)  of  24x24x16  cm3  for  an  individual  voxel 
volume  of  1.5  cm3;  256  bipolar  gradient  echo  pairs  (t  samples)  were  collected  with  bandwidth  of  1190  Hz.  A  non-water-suppressed  3D  EPSI  data  using  one  average 
was  used  for  eddy  current  and  phase  correction  of  the  suppressed  data.  With  pulse  repetition  time/echo  time  =  1500/30  ms  and  8  averages,  the  3D  EPSI  phantom  scan 
duration  was  12  minutes.  Global  water  suppression  was  performed  just  before  the  PRESS  localization. 

Results:  Both  the  non-water  suppressed  (NWS)  and  water  suppressed  (WS)  scans  were  separated  into  positive  (even)  and  negative  (odd)  subsets  first  and  reorganized 
into  kx-ky-kz-t  matrices.  Reconstruction  of  the  NUS  based  3D  EPSI  data  sets  was  done  offline  using  a  custom  MATLAB  software  package.  Fig.2.  shows  the  extracted 
ID  spectrum  from  the  central  voxel  from  the  prostate  phantom  (volume  of  1cm3).  We  were  able  to  detect  and  quantitate  glutamate/glutamine  (Glx),  spermine,  myo¬ 
inositol  in  addition  to  citrate,  choline  and  creatine.  The  metabolite  quantitation  was  done  by  Felix  peak  integral  quantitation  method.  The  mean  and  standard  deviation 
of  metabolites  with  respect  to  creatine  were:  Citrate  (12.642±3.230),  spermine  ((1.108±0.167),  choline  (0.619±0.069),  myo-inositol  (1.816±0.345)  and  Glx 
(2.239±0.267).  The  expected  and  calculated  experimental  metabolite  ratios  with  respect  to  creatine  are  shown  in  Fig.3. 


ADC 


Fig.l.  NUS  based  3D  EPSI  sequence  Fig.2.  An  extracted  spectrum  from  the  central  voxel 


Fig.3.  Comparison  of  experimental  and 
expected  metabolites  ratios  wrt  creatine 


Discussion:  By  using  CS  based  reconstruction,  we  have  demonstrated  a  potential  reduction  in  acquisition  time  by  up  to  80%  for  proton  MRSI,  with  negligible  loss  of 
information  as  evaluated  on  the  basis  of  clinically  relevant  metrics  .  The  NUS  3D  spatial  +1D  spectral  encoded  EPSI  sequence  recorded  short  TE-based  spectra  from 
multiple  regions  of  a  prostate  phantom  and  showed  more  metabolites  in  addition  to  the  normally  detected  citrate,  creatine,  choline  and  spermine. 

Conclusion:  The  NUS  based  3D  EPSI  data  using  CS  reconstruction  to  detect  and  quantify  metabolites  less  T2-weighted  than  the  earlier  3D  MRSI  sequences  using 
conventional  phase  encoding.  It  could  be  employed  to  help  detect  prostate  cancer  more  accurately  thereby  eradicating  negative  biopsies,  evaluate  cancer  stage 
noninvasively,  and  guide  treatment  selection.  Our  current  effort  is  focused  on  developing  an  optimal  sequence  with  interleaved  water-suppressed  metabolite  and 
unsuppressed  water  acquisition. 
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Abstract 

Overlap  of  metabolites  is  a  major  limitation  in  one-dimensional  (ID)  spectral-based 
single  voxel  MR  Spectroscopy  (MRS)  and  multivoxel-based  MR  spectroscopic  imaging 
(MRSI).  By  combining  echo-planar  spectroscopic  imaging  (EPSI)  with  two-dimensional 
(2D)  J  resolved  spectroscopic  sequence  (JPRESS),  2D  spectra  can  be  recorded  in 
multiple  locations  in  a  single  slice  of  prostate  using  four  dimensional  (4D)  Echo-Planar  J- 
Resolved  Spectroscopic  Imaging  (EP-JRESI)  combining  2  spectral  with  2  spatial 
dimensions.  The  goal  of  the  present  work  was  to  investigate  metabolite  levels  after  non¬ 
linear  reconstruction  of  an  accelerated  4D  EP-JRESI  data  using  Maximum  Entropy 
(MaxEnt)  and  Total  Variation  (TV)  in  Prostate  Cancer  (PCa).  Twenty  two  prostate  cancer 
patients  with  a  mean  age  of  63.8  years  (range:  46-79  years)  were  investigated  in  this 
study.  A  4D  non-uniformly  under-sampling  (NUS)  EP-JRESI  sequence  was  implemented 
on  the  Siemens  3T  MRI  scanner.  The  NUS  data  was  reconstructed  using  two  non-linear 
reconstruction  methods,  namely  MaxEnt  and  TV.  Significant  decrease  of  Cit  in  cancer 
locations  and  decreased  levels  of  Spermine  (Spm),  myo-Inositol  (ml),  and  increased 
levels  of  Glutamate  and  Glutamine  (Glx)  with  respect  (Choline  (Ch)+Creatine  (Cr))  were 
observed  in  cancer  locations  compared  to  non-cancerous  tissues  when  the  acquired  4D 
EP-JRESI  data  was  reconstructed  using  TV  and  MaxEnt.  We  have  shown  that  it  is 
possible  to  accelerate  the  4D  EP-JRESI  sequence  by  4X  and  the  data  can  be  reliably 
reconstructed  using  the  TV  and  MaxEnt  methods.  The  total  acquisition  duration  was  less 
than  13  minutes  and  we  were  able  to  detect  and  quantify  several  metabolites. 
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KEYWORDS:  Magnetic  Resonance  Spectroscopy,  Prostate  Cancer,  4D  EP-JRESI, 
Citrate,  Myo-inositol,  Glx,  Echo  planar  spectroscopic  imaging. 

Abbreviations  used:  PCa,  prostate  cancer;  PSA,  prostate  specific  antigen;  DRE,  digital 
rectal  examination;  MaxEnt,  maximum  entropy;  TV,  total  variation;  NUS,  non-uniform 
undersampling;  VOI,  volume  of  interest;  ROC,  receiver  operating  charecterstic;  PPV, 
positive  predictive  value;  NPV,  negative  predictive  value;  MRSI,  magnetic  resonance 
spectroscopic  imaging;  CS,  compressed  sensing;  JPRESS,  J-resolved  spectroscopy;  Ch, 
choline;  Cr,  creatine;  Cit,  citrate;  Spm,  spermine;  ml,  myo-inositol;  Glu,  glutamate;  Gin, 
glutamine;  ID,  one-dimensional;  2D,  two-dimensional;  3D,  three-dimensional;  4D,  four¬ 
dimensional;  RF,  radio-frequency;  SNR,  signal  to  noise  ratio;  EPSI,  echo  planar 
spectroscopic  imaging;  EP-JRESI,  echo-planar  J-resolved  spectroscopic  imaging; 
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Introduction 

Prostate  cancer  (PCa)  is  the  most  commonly  diagnosed  noncutaneous  malignancy  in  the 
USA  and  is  the  second-leading  cause  of  cancer-related  death  in  men  (1).  One  in  6  men 
will  be  diagnosed  with  PCa  during  his  lifetime  but  only  one  in  36  will  die  of  this  disease. 
Currently,  annual  prostate  specific  antigen  (PSA)  test  and  digital  rectal  examination 
(DRE)  are  routinely  performed  (2)  for  screening.  PSA  is  a  controversial  screening  test 
that  measures  the  serum  level  of  PSA  in  blood  samples,  as  65-75  %  of  PSA  screening 
gives  false-positive  results  leading  to  over-diagnosis  (3).  The  use  of  systematic  transrectal 
biopsy  can  miss  significant  prostate  cancers,  because  of  both  a  random  sampling  error  (4), 
and  because  up  to  a  third  of  significant  tumors  lie  in  the  anterior  part  of  the  gland,  based 
on  studies  of  radical  prostatectomy  specimens  (5).  Hence,  there  is  an  immediate  need  for 
early  yet  accurate  detection  of  PCa  to  improve  disease  outcomes. 

Proton  MR  Spectroscopy  (*H-MRS)  enables  detecting  a  range  of  biochemicals  in 
prostate  by  making  use  of  the  proton  signals  in  these  molecules.  Detection  of 
biochemicals  in  vivo  is  limited  to  concentrations  of  more  than  0.5-1  mM.  Signals  of 
citrate  (Cit),  creatine  (Cr),  choline  (Ch)  and  spermine  (Spm)  can  be  detected  throughout 
the  prostate,  with  increased  levels  of  Ch  and  decreased  levels  of  Cit  being  indicative  of 
cancer (6-8). 

Current  limitations  of  single  voxel  based  MRS  and  MR  spectroscopic  imaging 
(MRSI)  in  prostate  is  due  to  the  overlap  of  metabolite  resonances,  allowing  the 
quantification  of  only  a  few  metabolites  (Cit,  Ch,  Cr  and  Spm)  and  use  of  long  echo  times. 
The  conventional  MRSI  technique  can  be  accelerated  by  echo-planar  spectroscopic 
imaging  (EPSI)  (9-13).  The  EPSI  speeds  up  MRSI  by  using  an  echo-planar  read-out  of 
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one  spectral  and  one  spatial  dimensions,  thereby  achieving  an  acceleration  factor  equal  to 
the  number  of  points  along  one  of  spatial  dimensions.  For  example,  in  a  2D  spatial  matrix 
array  (16*16)  would  take  16  times  less  time  with  EPSI  than  the  conventional  MRSI. 
However,  the  acceleration  may  be  at  the  cost  of  signal  to  noise  ratio  (SNR)  (14)  and 
spectra  could  be  affected  by  Nyquist  ghost  artifacts  (15). 

A  single-voxel  based  2D  JPRESS  has  been  evaluated  in  PCa,  and  showed 
improved  spectral  dispersion  due  to  the  added  spectral  dimension  (16,17).  Compressed 
sensing  was  recently  introduced  as  a  powerful  method  to  accelerate  MRI  by  exploiting 
the  sparsity  of  the  images  in  a  known  transform  domain  to  reconstruct  non-uniformly 
undersampled  (NUS)  k-space  data  (18).  By  combining  EPSI  with  JPRESS,  2D  spectra 
can  be  recorded  in  multiple  locations  in  prostate  using  four  dimensional  (4D)  Echo- 
Planar  J-Resolved  Spectroscopic  Imaging  (EP-JRESI)  combining  2  spectral  with  2  spatial 
dimensions.  A  pilot  feasibility  was  demonstrated  recently  to  map  metabolites  in  the 
healthy  human  prostate  and  brain  (19,  20). 

The  Maximum  Entropy  method  (MaxEnt)  has  proved  to  be  an  enormously 
powerful  tool  for  reconstructing  images  from  many  types  of  MRI  data  but  it  has  been 
used  most  spectacularly  in  astronomy,  where  it  deals  routinely  with  images  of  up  to  a 
million  pixels,  with  high  dynamic  range  (21).  Maximum  entropy  and  total  variation  (TV) 
algorithms  have  been  used  to  reconstruct  NUS  indirect  spectral  and  spatial  dimensions 
(22). 

The  TV  algorithm  was  first  proposed  by  Rudin  et  al.  (23)  for  image  denoising  and 
since  then  has  been  successfully  used  for  image  restoration.  In  the  TV  algorithm,  an 
objective  function  using  the  TV  norm  is  minimized  subject  to  a  data  fidelity  posed  by  the 
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acquired  projection  data.  Minimizing  the  image  gradient  essentially  suppresses  those  high 
spatial  frequency  parts  such  as  streaking  artifacts  and  noise  in  the  reconstructed  images. 

The  goal  of  the  present  work  was  to  compare  MaxEnt  and  TV  reconstruction  of 
the  accelerated  4D  EP-JRESI  data  acquired  in  PCa  patients. 

Materials  and  Methods 

Patients:  Between  March  2012  and  May  2013,  twenty  two  PCa  patients  with  a  mean  age 
of  63.8  years  (range:  46-79  years)  who  underwent  radical  prostatectomy  were  selected 
for  the  study.  Patients’  Gleason  scores  varied  between  6  and  9.  Prostate-specific  antigen 
(PSA)  varied  from  0.7  to  22.8  ng/mL  (mean  of  6.23  ng/mL).  These  patients  were  scanned 
using  a  3T  Siemens  (Siemens  Medical  Solution,  Erlangen,  Germany)  MRI  scanner  with 
an  endorectal  "receive"  coil.  The  protocol  combining  MRI  and  MRS  was  performed  at 
least  8  weeks  after  transrectal  ultrasound-guided  sextant  biopsy.  The  entire  protocol  was 
approved  by  the  Institutional  Review  Board,  and  informed  consent  was  obtained  from 
each  patient.  Prostate  cancer  was  histopathologically  confirmed  and  graded  after  radical 
prostatectomy.  The  voxels  covering  the  tumorous  lesions  from  the  peripheral  zone  (PZ) 
were  selected  and  indicated  as  tumor  voxels.  After  reconstruction  the  EP-JRESI  were 
overlaid  onto  MRI  images. 

MRI  and  MRSI:  A  combined  body  matrix  phased-array  coil  assembly  and  an  endorectal 
coil  were  used  in  the  "receive"  mode  while  a  quadrature  body  coil  was  used  for 
"transmit".  All  patients  were  imaged  in  supine  feet- first  position.  Axial  images  were 
oriented  to  be  perpendicular  to  the  long  axis  of  the  prostate,  which  was  guided  by  the 
sagittal  images.  Axial,  coronal,  and  sagittal  T2-weighted  (T2W)  turbo  spin-echo  imaging 
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was  performed  with  the  following  parameters:  repetition  time/echo  time  (TR/TE),  3850- 
4200/96-101  ms;  section  thickness,  3  mm;  field  of  view,  20  x  20  cm  ;  echo  train  length, 
13;  matrix,  320  x  256. 

A  4D  NUS  EP-JRESI  sequence  was  implemented  on  the  3T  MRI  scanner  and  the 
volume  of  interest  (VOI)  was  localized  using  three  slice-selective  radio-frequency  (RF) 
pulses  (90°-180°-1800)  (Fig.l).  The  total  time  for  obtaining  a  fully  sampled  4D  EP-JRESI 
scan  (TR  of  1.5s,  16ky*16kx,  64-100ti,  512t2)  can  be  more  than  25  minutes.  The 
parameters  for  the  EP-JRESI  were:  TR/TE/Avg  =  15 00/3 0ms/ 1-2,  16  phase  encoding 
steps,  512  complex  points  with  an  F2  bandwidth  of  1190Hz.  For  the  second  dimension 
(Fi),  64  increments  with  bandwidths  of  1000Hz  were  used.  A  4X  NUS  was  imposed 
along  the  plane  containing  incremented  spectral  and  spatial  dimensions  (ti  and  ky).  The 
individual  voxel  volume  in  human  prostate  was  1ml.  Two  sets  of  data  were  collected,  one 
water  suppressed  scan  (WS)  with  a  total  scan  time  of  6-12  minutes  and  a  second  non¬ 
water  suppressed  scan  (NWS)  using  one  average  only. 

Data  Analysis:  The  NUS  data  were  reconstructed  by  MaxEnt  and  TV  separately.  A 
modified  Split-Bregman  algorithm  (21)  solves  the  unconstrained  TV  optimization 
problem  as  shown  in  Eqn.(l)  below: 

min  ||  V m||i  +k||Fum-y||2  (1) 

m 

where  V  is  the  gradient  operator,  m  is  the  reconstructed  data,  ||x||i  is  the  1 1  -norm,  X  is  a 
regularization  parameter,  Fu  is  the  under  sampled  Fourier  transform,  and  y  is  the  under¬ 
sampled  data.  The  above  equation  removes  the  incoherent  artifacts  due  to  NUS  by 
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minimizing  the  total  variation  (TV)  while  maintaining  consistency  with  the  sampled 
measurements. 

MaxEnt  is  a  constrained  convex  optimization  algorithm  that  uses  a  variant  of  the 
conjugate  gradient  method  to  iteratively  solve  the  inverse  problem  (21,  24,  25): 

maximize  S1/2  (f)  s.t.  HF'1  Kf-D||2  <  a  (2) 

where  /  is  the  estimated  fully-sampled  spectrum  at  each  iteration,  F1  is  the  inverse 
Fourier  transform,  K  is  the  NUS  matrix,  D  is  the  measured  time-domain  data,  a  is  the 
noise  standard  deviation,  and  S  1/2(0  is  the  spin- 14  entropy  of  the  estimated  spectrum  (26). 
All  NUS  based  4D  EP-JRESI  data  were  processed  using  TV  and  MaxEnt  reconstruction 
with  custom  MATLAB  software. 

Statistical  Analysis: 

Statistical  analyses  were  performed  with  SPSS  21  (SPSS  Inc.,  Chicago,  IL,  USA).  Using 
logistic  regression  analysis,  areas  under  the  curve  (AUC)  of  the  receiver  operating 
characteristic  (ROC)  were  calculated  for  various  metabolites  in  discriminating  between 
MaxEnt  and  TV  reconstruction  methods.  Also  the  paired  t-test  was  used  to  determine  the 
various  metabolite  ratios  in  cancer  and  non-cancer  locations.  A  p-value  less  than  0.05 
was  considered  statistically  significant. 

Results 

We  were  able  to  detect  Cit,  Ch,  Cr,  Spm,  myo-inositol  (ml),  and  glutamate  plus 
glutamine  (Glx)  in  cancer  and  non-cancer  locations  and  quantified  using  the  peak 
integration  method.  Figs.  2  and  3  show  the  metabolite  ratios  of  Cit,  Spm,  ml,  and  Glx  of 
cancer  and  non-cancer  locations  processed  by  MaxEnt  and  TV.  The  mean  metabolite 
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ratios  and  standard  deviations  (SD)  of  Cit,  Spm,  ml  and  Glx  with  respect  to  (Ch+Cr) 
were  quantified  for  the  TV  reconstructed  non-cancer  locations  as:  1.215  ±0.537,  0.801 
±0.622,  0.214  ±0.106,  0.228  ±0.105.  In  the  cancer  location,  the  metabolites  ratios  were: 
0.361  ±0.241,  0.607  ±0.457,  0.202  ±0.073,  0.246  ±0.146.  Similarly  the  mean  metabolite 
ratios  and  SD  of  Cit,  Spm,  ml  and  Glx  with  respect  to  (Ch±Cr)  were  quantified  for  the 
MaxEnt  reconstructed  non-cancer  locations  as:  1.234  ±0.557,  0.929  ±0.629,  0.238  ±0.101, 
0.219  ±0.105.  In  the  cancer  locations,  the  metabolite  ratios  were:  0.347  ±0.188,  0.581 
±0.335,  0.210  ±0.144,  0.234  ±0.175. 

We  found  that  mean  Cit  metabolite  ratios  were  significantly  lower  in  cancer  locations 
compared  to  non-cancerous  locations  in  both  TV  and  MaxEnt  reconstruction.  Decreased 
levels  of  Spm,  ml  and  increased  levels  of  Glx  with  respect  to  (Ch±Cr)  ratio  were 
observed  in  cancer  locations  compared  to  non-cancer  locations  in  the  TV  and  MaxEnt 
reconstructions.  There  were  no  statistically  significant  changes  when  metabolites  ratio 
were  correlated  with  GS. 

Fig.4  shows  a  spatial  maps  of  (Ch±Cr)  for  the  TV  and  MaxEnt  reconstructed  data 
acquired  in  a  74  year  old  PCa  patient.  The  MaxEnt  and  TV  reconstructions  of  cancer 
(B&D)  and  non-cancer  (C&E)  locations  extracted  from  the  JPRESS  spectrum  of  a  4D 
EP-JRESI  voxel  is  shown  in  Fig.5.  We  have  compared  and  correlated  the  TV  with 
MaxEnt  reconstruction  methods  of  Cit,  Spm,  ml  and  Glx  in  the  cancer  and  non-cancer 
locations.  The  correlation  of  the  Cit/(Ch±Cr),  Spm/(Ch±Cr),  mI/(Ch±Cr)  and 
Glx/(Ch±Cr)  ratios  for  the  MaxEnt  and  TV  reconstructions  in  cancer  and  non-cancer 
locations  are  shown  in  Fig.6.  A  positive  correlation  was  found  in  the  following 
metabolites  in  cancer  locations:  Cit/(Ch±Cr)  (R2=0.963),  Glx/(Ch±Cr)  (R2=0.918), 
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Spm/(Ch+Cr)  (R2=0.905)  and  mI/(Ch+Cr)  (R2=0.820).  Also  we  observed  positive 
correlation  of  the  above  mentioned  metabolites  in  non-cancer  location  (not  reported). 

Detailed  results  of  the  logistic  regression  analysis  and  consequent  receiver 
operating  characteristic  (ROC)  curve  analyses  are  given  in  Table  1  including  the 
sensitivity,  specificity,  positive  predictive  value  (PPV),  negative  predictive  value  (NPV), 
area  under  curve  (AUC),  and  accuracy  for  classifying  the  MaxEnt  and  TV  methods.  ROC 
curve  analyses  for  differentiating  the  metabolite  ratios  of  cancer  and  non-cancer  locations 
in  MaxEnt  suggest  Cit  gives  the  best  predictability  with  sensitivity  86.4%,  specificity 
86.4%,  accuracy  86.4%  and  AUC  94.0%.  Also  the  overall  sensitivity,  accuracy  and  AUC 
of  ml  and  Spm  were  slightly  higher  in  MaxEnt  compared  to  TV. 

Discussion 

Using  NUS  with  non-linear  iterative  reconstruction,  we  have  demonstrated  a  potential 
reduction  of  acquisition  times  by  up  to  75%  for  EP-JRESI  with  minimal  loss  of 
information.  In  addition  to  significant  decrease  of  Cit  in  cancer  locations,  decreased 
metabolite  ratios  of  Spm,  ml  and  increased  Glx  were  found  in  cancer  locations  compared 
to  non-cancer  locations.  Even  though  TV  and  MaxEnt  reconstruction  showed  comparable 
results  in  cancer  and  non-cancer  locations,  there  was  a  slight  increase  of  sensitivity, 
accuracy  and  AUC  in  MaxEnt  reconstruction. 

In  the  present  study,  we  report  the  relative  concentrations  of  Cit,  Spm,  ml  and  Glx 
with  respect  to  (Ch+Cr),  because  the  total  Cr  peak  (3.0  ppm)  is  very  close  to  the  total  Ch 
peak  (3.2  ppm)  in  these  in  vivo  MR  spectra  and  may  not  always  be  separable.  In  this 
study,  significantly  higher  level  of  Cit  was  observed  in  the  peripheral  zone  of  healthy 
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prostates.  The  production  and  storage  of  citrate  is  one  of  the  main  characteristics  of 
normal  prostate. 

It  is  likely  that  the  drop  in  Cit  precedes  malignant  transformation  (27).  It  has  been 
suggested  that,  due  to  a  metabolic  switch,  neoplastic  cells  would  become  Cit  oxidizing, 
unlike  normal  prostatic  cells  that  show  a  low  Cit-oxidizing  capability  (28).  Decreased 
levels  of  zinc,  which  would  relieve  m-aconitase  from  inhibition,  has  been  proposed  as 
one  of  the  reasons  for  the  decreased  level  of  Cit  in  PCa  (27). 

Choline  is  an  essential  component  of  cell  membrane  synthesis  and  phospholipid 
metabolism  and  functions  as  an  important  methyl  donor.  Choline-containing  molecules 
are  an  essential  component  of  cell  membranes,  which  are  more  highly  concentrated  in 
tumorous  areas  within  the  prostate  than  in  healthy  prostate  tissue.  (29,30).  Choline  groups 
are  precursors  and  breakdown  products  of  the  phospholipid  phosphatidylcholine,  a  major 
cell  membrane  compound  (31).  Increased  Ch  is  observed  due  to  an  altered  phospholipid 
metabolism  in  PCa  cell  lines  (30).  This  alteration  is  most  likely  due  to  an  increased 
expression  and  activity  of  choline-kinase,  a  higher  rate  of  choline  transport,  and  an 
increased  phospholipase  activity  (31). 

The  polyamines  spermine,  spermidine,  and  putrescine  are  essential  for  the 
differentiation  and  proliferation  of  cells,  the  synthesis  of  DNA,  RNA,  and  proteins,  and 
the  stabilization  of  cell  membranes  and  cytoskeletal  structures  (32).  Previous  studies 
observed  high  levels  of  spermine  in  healthy  prostate  tissue  and  BPH  and  reduced 
spermine  levels  in  malignant  prostate  tissue  (17,  33-35). 

The  osmoregulator  myo-inositol  is  expressed  in  a  variety  of  tissues,  and  its 
decrease  was  observed  in  PCa  within  human  expressed  prostatic  secretions  (EPS)  using 
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high  resolution  NMR  (36)  and  in  breast  tumors  (37).  In  our  study,  slightly  decreased  ml 
were  observed  in  cancer  locations  but  were  not  statistically  significant. 

Glu  and  Gin  are  difficult  to  be  resolved,  owing  to  the  overlap  of  resonances.  As  a 
result,  most  reliable  MRS  studies  use  the  sum  of  Glu  and  Gin  (expressed  as  either  Glx  or 
Glu  +  Gin).  Glutamate  is  extensively  involved  in  metabolic  and  oncogenic  pathways. 
Koochekpour  (38)  showed  that  serum  glutamate  levels  directly  correlated  with  Gleason 
scores  (<6  vs.  >8)  and  primary  PCa  aggressiveness.  In  our  study,  elevations  of  Glx  were 
found  in  the  cancer  locations  but  were  not  statistically  significant. 

There  was  an  overlap  between  cancer  and  non-cancer  locations,  possibly  due  to 
the  low  signal-to-noise  ratio  of  the  Cr  peak.  In  addition,  Ch  and  Cr  were  not  always  well 
resolved  especially  in  cancer  locations,  adding  to  the  uncertainty  in  quantitation. 
However,  using  prior  knowledge  fitting  (ProFit),  accurate  quantitation  of  metabolites  can 
be  improved  (39).  Due  to  the  small  patient  population,  we  did  not  find  any  significant 
changes  in  other  metabolites.  Also  this  study  focused  in  the  peripheral  zone  of  the 
prostate  where  80%  of  cancer  occurs.  The  advantage  of  NUS  based  4D  EP-JRESI 
sequence  is  in  recording  short  TE-based  spectra  from  multiple  regions  of  human  prostate 
and  more  metabolites,  such  as  myo-inositol,  spermine  and  Glx  in  addition  to  the  normally 
detected  citrate,  creatine  and  choline. 

Conclusion 

We  were  able  to  detect  metabolites  using  NUS  based  4D  EP-JRESI  data  acquired  in 
clinically  acceptable  times  (<12  minutes).  We  have  shown  that  it  possible  to  undersample 
the  4D  EP-JRESI  sequence  with  an  acceleration  factor  of  4X  and  the  data  can  be  reliably 
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reconstructed  using  the  TV  and  MaxEnt  methods.  Both  non-linear  reconstruction 
methods  provided  comparable  results. 
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Table  1.  Measures  of  Sensitivity,  Specificity,  PPV,  NPV,  and  Accuracy  of  MaxEnt 
and  TV  Method  using  ROC  Curve  Analysis 


MaxEnt 

Metabolites 

Sensitivity  % 

Specificity  % 

Accuracy  % 

PPV  % 

NPV  % 

AUC  % 

Cit 

86.4 

86.4 

86.4 

86.4 

86.4 

94.0 

Glx 

40.9 

59.1 

50.0 

50.0 

50.0 

48.0 

ml 

63.6 

63.6 

63.6 

63.6 

63.6 

65.0 

Spin 

81.8 

40.9 

61.4 

56.7 

64.3 

64.9 

TV 

Cit 

86.4 

86.4 

86.4 

86.4 

86.4 

92.1 

Glx 

54.5 

47.6 

51.2 

52.2 

50.0 

50.6 

ml 

50.0 

38.1 

44.2 

45.8 

42.1 

50.6 

Spin 

77.3 

40.9 

59.1 

58.1 

69.2 

58.5 
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Figure  Legends 

1.  Fig.l.  Schematic  diagram  of  4D  EP-JRESI  pulse  sequence. 

2.  Fig.2.  Metabolite  ratios  of  Cit,  Spm,  ml  and  Glx  in  cancer  and  non-cancer  locations 
processed  by  MaxEnt. 

3.  Fig.3.  Metabolite  ratios  of  Cit,  Spm,  ml  and  Glx  in  cancer  and  non-cancer  locations 
processed  by  TV. 

4.  Fig.4.  Spatial  map  of  (Ch  +Cr)  for  TV  and  MaxEnt  reconstructions  of  4D  EP-JRESI 
data  recorded  in  a  74  year  old  PCa  patient. 

5.  Fig.5.  2D  JPRESS  spectra  extracted  from  an  EP-JRESI  voxel  (A)  after  the  MaxEnt 
and  TV  reconstructions  of  cancer  (B&D)  and  non-cancer  (C&E)  locations. 

6.  Fig.6.  Correlation  of  MaxEnt  and  TV  for  Cit,  Glx,  ml  and  Spm  in  cancer  locations. 
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Fig.l.  Schematic  diagram  of  4D  EP-JRESI  pulse  sequence. 
116x58mm  (300  x  300  DPI) 
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Fig. 2.  Metabolite  ratios  of  Cit,  Spm,  ml  and  Glx  in  cancer  and  non-cancer  locations  processed  by  MaxEnt. 
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Fig. 3.  Metabolite  ratios  of  Cit,  Spm,  ml  and  Glx  in  cancer  and  non-cancer  locations  processed  by  TV. 

157x115mm  (300  x  300  DPI) 
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TV  MaxEnt 

Fig. 4.  Spatial  map  of  (Ch  +Cr)  for  TV  and  MaxEnt  reconstructions  of  4D  EP-JRESI  data  recorded  in  a  74  year 

old  PCa  patient. 
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Fig. 5.  2D  JPRESS  spectra  extracted  from  an  EP-JRESI  voxel  (A)  after  the  MaxEnt  and  TV  reconstructions  of 

cancer  (B&D)  and  non-cancer  (C&E)  locations. 

148x126mm  (300  x  300  DPI) 
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